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Abstract. We make a direct and model-independent measurement of the low tt+tt^ mass phase 
motion in the D+ n^n^n^ decay. Our preliminary results show a strong phase variation, 
compatible with the isoscalar (7(500) meson. This result confirms our previous result [1] where 
we found evidence for the existence of this scalar particle using full Dalitz-plot analysis. We apply 
the Amplitude Difference (AD) method [2] to the same Fermilab E791 data sample used in the 
preceding analysis. We also give an example of how we extract the phase motion of the scalar 
amplitude, looking at the /o(980) in — > n^n'^n'^ decay. 

INTRODUCTION 

Fermilab experiment E791, with a full Dalitz plot analysis, showed strong evidence for 
the existence of light and broad scalar resonances in charm meson decay [1,3]. The 
n^n^ resonance is compatible with the isoscalar meson a(500), and was observed in 
the Cabbibo-suppressed decay k^k^tz^ . To get a good fit quality in this analysis, 

it was necessary to include an extra scalar particle, other than the well established dipion 
resonances [4] . For the new state, modeled by a Breit-Wigner amplitude, it was measured 
a mass and width of 478^23 ± ^'^ MeV/c^ and 324lJo ± 21 MeV/c^respectively . The 
— i> a(500);r+ decay contribution is dominant, accounting for approximately half 
of this particular D+ K^n^n^ decay. We found also evidence for a scalar K~n^ 
resonance, or K, in the Cabibbo-allowed decay K^n^n^ [3]. Further studies 

about K are discussed in this proceeding [5]. 

In full Dalitz plot analyses, each possible resonance amplitude is represented by a 
Breit-Wigner function multiplied by angular distributions associated with the spin of 
the resonance. The various contributions are combined in a coherent sum with complex 
coefficients that are extracted from maximum likelihood fits to the data. The absolute 
value of the coefficients are related to the relative fraction of each contribution and the 
phases take into account the final state interaction (FSI) between the resonance and the 
third pion. 

Due to the importance of this scalar meson in many areas of particle and nuclear 
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physics, it is desirable to be able to confirm the amplitude's phase motion in a direct 
observation, without having to assume, a priori, the Breit-Wigner phase approximation 
for low-mass and broad resonances [6, 7, 8]. Recently, a method was proposed to 
extract the phase motion of a complex amplitude in three body heavy meson decays 
[2]. The phase variation of a complex amplitude can be directly revealed through the 
interference in the Dalitz-plot region where it crosses with a well established resonant 
state, represented by a Breit-Wigner. 

Here we begin with a simple example, showing that the AD method can be applied to 
extract the resonant phase motion of the scalar amplitude due to the resonance /o(980), 
using the same /o(980) resonance in the crossing channel in the Dalitz plot of the decay 
Df — > ;r~7r"'"7r"'" using E791 data [9]. This example shows the ability of this method to 
extract the phase motion of an amplitude. Then we apply the AD method using the well 
known /2(1270) tensor meson in the crossing channel, as the base resonance, to extract 
the phase motion of the scalar low mass nn amplitude in D+ n~n~^n~^ , confirming 
the a (500) suggested by the E791 full Dalitz plot analysis [1]. 

EXTRACTING /o(980) PHASE MOTION WITH THE AD 

METHOD. 

From the original 2 x 10^^ event data collected in 1991/92 by Fermilab experiment E791 

from SOOGeV /c K~ —nucleon interactions [10], and after reconstruction and selection 
criteria, we obtained the n'K^K^ sample shown in Figure 1. To study the resonant 
structure of these three-body decays we consider the 1686 events with invariant mass 
between 1.85 and 1.89 GeV/c^, for the Z)+ analysis and the 937 events with invariant 
mass between 1.95 and 1.99 GeV/c^ for the D+. Figure 2(a) shows the Dalitz-plot for 
the Df 7U~7U~^7t~^ selected events and Figure 2(b) the Dalitz-plot for Z)+ 7C~7U~^7t~^ 
events. The two axes are the squared invariant-mass combinations for 7U~7C~^, and the 
plot is symmetrical with respect to the two identical pions. 

We can see in Figure 2(a) the scalar /o(980) in 512, the square invariant mass, crossing 
the /o(980) in 513, forming an interference region around 513 = 512 = 0.95GeV^. The AD 
method uses the interference region, between two crossing resonances, to extract the 
phase motion of one of them, and Final State Interaction (FSI) phase, provided that the 
second is represented by a Breit-Wigner [2]. In fact we are using a Bootstrap approach; 
that is, using a well established resonance /o(980) in 512 to extract its phase motion in 
513. It is a nice and didactic example to show the ability of this method to extract the 
phase motion of an amplitude and the FSI phase within the E791 data sample. 

The coherent amplitude to describe the crossing between a well known scalar reso- 
nance, represented by a Breit-Wigner in sn, and a complex amplitude under study in 513 
in a limited region of the phase space, where we can neglect any other contributions, is 
given by: 

■c/(^12,*13) = ai?^^(*12)+a./(;?7v^l3)*'«5(5l3)e'(^(-^'^)+^) (1) 

P* /V^i3 is ^ phase space factor to make this description compatible with nn scatter- 
ing, 7 is the final state interaction (FSI) phase difference between the two amplitudes. 
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FIGURE 1. The n n^n^ invariant mass spectrum. The dashed line represents the total background. 
Events used for the Dalitz analyses are in the hatched areas. 
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FIGURE 2. (a) The D+ -> 7r"7r+7r+ DaHtz plot and (b) the D+ n'n+n+ Dahtz plot. Since there 
are two identical pions, the plots are symmetrical. 



qr and tts are respectively the real magnitudes of the resonance and the under-study 
complex amplitude. Finally sin5{s\^)e^^'^'^^^^ represents the most general amplitude for a 
two-body hadronic interaction. 

The Breit Wigner distribution is given by: 




FIGURE 3. /o(980) sn distribution, divided in mo + e (red) and mo — e (blue) 
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Taking the amplitude square of Equation 1 we get: 



*12,*13) 1^=4 I =^^/o(980)l 



(512) 1^ +a^/p*^/si3 sin^d{si3 
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13) 



{ml -si2)^ + mlr^{s 12) 



X [{ml-si2)cos{5{si3)+r)+morosin{5{si3)+r)] (2) 



Since the Breit-Wigner is approximately symmetrical around niQ as seen in Figure 3 
(the asymmetries would come from r(m), and is negligible for the narrow /o(980)). We 
can divide our /o(980) mass distribution in two pieces, one for mq + £ and the other with 
niQ — £. From Equation 2 and noticing that the pure Breit-Wigner term will cancel we 
can write: 



Only the real part of the interference term in Equation 2 remains. 

To extract the phase motion of the scalar amplitude in 513 through the /o(980) in sn, 
represented by a Breit-Wigner, we took the events in sn between 0.7 and 1.2 GeV^ and 
divided them into two bins, as presented in Figure 3. The S13 distribution for the events 
of the region integrated between 0.95 and 1 .2 GeV^, is shown in Fig. 4a and the same 
in Figure 4b for events integrated between 0.7 and 0.95 GeV^. 
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FIGURE 4. si3 distribution, a) For events /'"" ^ \ ^(suj^b) P dsi2. b) For events 
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We can see that the peaks in these two plots are in different positions. The 
subtraction of these distributions, corresponds to the integration of Equation 3, that we 
can write as: 





2^(^12,^13) \^ dsn- 
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ds 
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^(5m(25(5i3) + 7) -siny) (4) 

where is a constant factor coming from the constant and integrated factors of Equation 
3, to be determined from data. The variation of the phase space in the integral was 
considered negligible for the /o(980) resonance. A / s^^ directly reflects the behavior 
of 5(5'i3). A constant A | ,2/ p would imply constant 5(513). This would be the case for 
a non-resonant contribution. The same way a slow phase motion will produce a slowly 
varying A | ^ p and a full resonance phase motion produces a clear signature in A | p 
with the presence of zero, maximum and minimum values. 

The subtracted distribution, corresponding to Equation 4, is shown in Figure 5. There 
is a significant difference between the minimum (bin3) and maximum (bin4) of A / 

We can see in Equation 4 that the zeros occur when 5{sn) = 0*^, 180^ or ;r/2 — 7. In 
Figure 5 we can see a zero at 513 near 0.5GeV^, another one at = 0.95GeV^ and a 
third zero near 1.4GeV^. Assuming 5(513) is an analytical function of 513, Equation 4 
allow us to set the two following conditions at the maximum and minimum values of 
A / respectively: 

A / .^L - ^'■"(25(513) + 7) = -1 (5) 
A/<L^^m(25(5i3) + 7) = l (6) 



FIGURE 5. 513 distribution of A/ si^'^dsn. 



with these two conditions we get ^ and 7, calculated from the maximum and mini- 
mum values of the A / distribution in Figure 5: 

^=(a|^^-a|<2j/2 (7) 

.„-i.A^<,+ A/^ 
^ ^A/<.„-A/<2^J 

From Figure 5 and using the equations above, we measure 7= —0.15 ±0.31, that is 
compatible with zero, as should be since we are crossing the same resonances with, of 
course, the same final state interaction phase. 

With the above conditions we solve Equation 4 for 5(si3): 

5(513) = \{sin-\]^^ I ^(513) |2 +5m(7)) - 7) (9) 

Assuming that 5(^13) is an increasing function of 513, we can extract directly the 
5(513) value from each bin of Figure 5, creating the /o(980) phase motion shown in Fig- 
ure 6. The errors in the plot were produced by generating statistically compatible exper- 
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iments, allowing each bin of \ ^(512,513) p (Figure4a) and ^ \ .2/(512,513) p 

(Figure4b) to fluctuate randomly following a Poisson law. We then solve the problem for 
each "experiment". The error in each bin of 5 (513) will be the RMS of the distributions 
generated by the "experiments". 

From Figure 6 we can see what one could expect, that is the scalar amplitude near 
970GeV with a phase motion of about 180° degrees. This example demonstrates the 
ability of AD method to extract the phase motion of an amplitude with E791 statistics. 
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FIGURE 6. 5(ii3) plot with the errors. 

EXTRACTING THE SCALAR LOW MASS tztz AMPLITUDE 
PHASE MOTION WITH THE AD METHOD. 

In the preceding section, we showed how to apply the AD method to extract the phase 
motion of an amplitude, from nonleptonic charm-meson three-body decay. Here we 
apply the same method to extract the phase motion of the scalar low-mass mt amplitude 
in D+ — > %~Tt'^Tt^ decay, where we previously found strong experimental evidence 
for the existence of a light and broad isoscalar resonance [1]. To start this analysis, we 
have to decide what is the best well-known resonance to be used for crossing the low 
mass amplitude under study. Taking a look at Figure 2b we can see the signature of 
three resonances that in principle could be used, the p(770), /o(980) and /2(1270). In 
fact, the E791 analysis of this Dalitz plot found a significant contribution from these 
three resonances in n~7t~^n^ decay [1]. Since this D"*" decay is symmetric for 

the exchange of the tC^ meson, each resonance in su is present also in 513. Then if we 
use p(770) as the base resonance in 512, we have also the presence of the p(770) in 
same mass square distribution of the cj(500) in 513. The proximity of the p(770) with 
the a (500), both broad resonances, creates an overlap between them such that we are 
not able to separate the phase motion of one from the other. We could use the /o(980) 
as a base resonance, but again the presence of the p(770) overlapping with the <7(500) 
creates the same problem. 

There remains only the tensor meson /2(1270) candidate at Mq = l.6lGeV^, which is 
placed where the p (770) contribution reaches a minimum due the angular distribution 
in the middle of the Dalitz plot, as we can see from the Z)+ — > p(770);r+ decay Monte 
Carlo simulation shown in Figure 7. 

The amplitude for the crossing of the /2(1270) in S12 and the complex amplitude 
under study in i' 13 is given in the same way as in Equation 1: 




FIGURE 7. MC p (770) distribution in D+ n+n+ decay. There is little contribution between 1 .2 
to 1.8GeV2 



^(^12,^13) =6!^; <^^/2(1270)('^12) ^ ^^^2(1270) (■^12, ^13) + (10) 

+ a,/(p7v/^i3) sin5{su) e'(^(-^>^)+^) 

where ■'^^-^/2(i270)('^i27>^i3) is the angular function for the /2(1270) tensor resonance. 
The amplitude under study represents the scalar low mass kk amplitude in a limited 
region of the phase space, where we can neglect the other amplitude contributions. 

Both the width r{si2) and the angular function ■'^^^/2(i270) from this resonance 
produce asymmetries in 512 and consequently we can not use the nominal /2(1270) 
mass to divide our sample into two slices, as we did for the /o(980) example. So we 
performed a Monte Carlo study to determine the effective mass we must use. The S12 
Monte Carlo projection of the /2(1270) in D+ — > n^n^n^ decay is shown in Fig. 8. We 
can see the asymmetry created around the nominal /2(1270) mass value due to r(si2) 
and •'^^-^/2(i270) contributions to the amplitude. 

Here we require an effective mass squared (nieff), such that the number of events 
integrated between fn^j^j: and m^ff + £ is equal, by construction, to the number of events 
integrated between mg^y and m^yy — £. We choose, using the /2(1270) Monte Carlo 
distribution, a mass of w^yy = l.535GeV^, within ±0.26GeV^ ^, in such way that we 
can write: 



Within this mass region, the amount of p (770) events was estimate to be around 5% 
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/ , I ^^/2(1270) (^12)^=^^/2(1270) l^dsn = I ^^/2(1270)(^12)^'=^^/2(1270) 

(11) 

The effective mass squared meff and the separation between m^^j + £ (red) and 
m^yy — £ (blue) are shown in Figure 8. 

The ■'^^^/2(i270) function in sij, is presented in Figure 9^. The distribution between 
m^eff ^'^^ '"e// + £ is shown in Fig. 9a, for events between m^yy and m^yy — £ in Fig. 
9b. We can see that these two plots are slightly different. However we considered the 
approximation ■^'"^'^y^(i270) (■^13-' ~ ''^^■^/2(1270) (■^13-' ^'^^ average function 

^^^•y£f2[\2'jQ){si-i). Another important effect, that we had to take into account, is the 

zero of this function at = OASGeV^. Below we discuss the consequences of that in 
the AD method. 

With the above considerations about the /2(1270) in si2 and 513 we can write the 
integrated amplitude-square difference as: 



^ dsu 
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-^isin{2d{si3) + r)-sinr) ^'"^^/2(1270) (^13)/(p7v^13) (12) 



^ Since we divided our data sample by this function, we represent this function in a histogram with the 
same binning of data. 
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FIGURE 9. Fast MC '"^-#/2(i270) distribution in ^13. a) For events between m^^^ and in^j-j- + e. b) For 
events between m^^^ and rn^j-j- — e. 

This Equation is similar to Equation 4, with an extra angular function term 

The background and the acceptance are similar between m^ff and m^ff + £ and «^gyy 
and mgyy — £. Since we are subtracting these two distributions, we do not take into 
account these effects in our analysis. 

The / =2/2 in 513, for events integrated in 512 between w^yy = 1 .535GeV^ and w^yy + e 

and m^yy and w^yy — e, with e = 0.26Gey2 are presented in Figure 10a and b respec- 
tively. 

Subtracting these two histograms, in the same way we did for the /o(980) example, 
gives the A / of the Equation 12. The result is shown in Figure 11. 

Here we can not extract directly the phase motion from Figure 1 1 , as we did for the 
/o(980) example using the conditions 5 and 6. We have to divide the A / by 
(average of the distributions in Figures 9a and b) and multiply by p', since phase space 
here is an important effect. By doing this the only 513 dependence of the right hand 
side of Equation 12 is in 5(513). However, as we could see in Figure 9, there is a zero 
about 0.48GeV2 in the angular function, which means a singularity around this value 
in A/ To avoid this singularity, we first produced a binning in such a way that 
the singularity is placed in the middle of one bin. In Figure 12 we show the A / =2/^ by 
^ distribution. We can see that the 6th bin (around 0.48GeV2), has a huge error, that 
corresponds to the bin of the singularity. Due to the singularity we decided not to use 



^ For short we shall use, from here on ^.^2(1270) (-^is) = P*/V^i3 = P'- 




FIGURE 10. Events distributions in si^, a) for events f'"" ^ \ ^(ini^is) P dsn, b) For events 

/;|_J^(.12,.13)N^12. 



20 - 



10 - 



FIGURE 11. ii3 disti-ibution for A / 



this region (bin 6) further in this analysis. The consequences of this choice are going to 
be taken care of in systematic error studies. In any case, the singularity could only affect 
the position of the minimum of Figure 12. It does not interfere with the general feature 
of starting at zero, having statistically significant maximum and minimum values, and 
coming back to zero, indicating a strong phase variation. Bins 2 and 5 are respectively 
the maximum and minimum value of A / £/^p' of Figure 12 where we use the 
Equation 5 and 6 conditions. 

With the same assumptions used for the /o(980), that is 5(513) is an analytical and 
increasing function of 513, and using Equation 7, 8 and 9 (multiplied by p' and divided 
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FIGURE 12. ii3 distribution for A / pi 

by we can extract y and 5(513) from Figure 12. For the FSI phase we found 
7 = 3. 26 ±0.33, that is somewhat bigger than found by the E791 full Dalitz-plot analysis 
[1] {yDaiitz = 2-59 ± 0.19). The fact that we used the effective mass for the /2(1270) = 
1.535 GeV^ instead of the nominal mass is responsable for the shift observed in the 
relative phase. To verify this statement we generated 1000 samples of fast-MC, with 
only two amplitudes, /2(1270) and <7(500). For both we used Briet-Wigner functions 
and the E791 parameters. We generated the phase difference of 2.59 rad, measured by 
the E791. For these 1000 samples, we measure 7 using the method presented here. The 
result has a mean value of 3.07. We can say that the difference between the generated 
and measured 7 value is a correction factor due to the use of an effective mass. Using 
this offset factor ( 2.59 - 3.07 = -0.48) we correct the measurement 7 ~ 3.26 ± 0.33 to 
Ycorr = 2.78 ±0.33. So the observed 7 difference between Dalitz analysis and the Ycorr 
is in good agreement, with a difference below one standard deviation. 

The 5(513) was extracted bin by bin, with the same approach for the errors used in the 
/o(980) example, and we got the distribution shown in Figure 13 ^. We can see a strong 
phase variation of about 180^ around the mass for the (7(500), showing a phase motion 
compatible with a resonance. 



CONCLUSIONS 

We showed that the AD method can be applied to E791 data to extract the phase motion 
of the resonance /o(980) in the Dalitz plot of the decay Df — > 7t~ Tf^Tf^ . This example 
demonstrates the ability of this method to extract the phase motion of a resonance 



In Figure 13 there is no the 6th bin because of the singularity we mentioned above. 
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FIGURE 13. Phase motion 5(si3) distribution for the scalar low mass nn amplitude, with the errors. 



amplitude. 

Preliminary E791 results present a direct and model-independent approach, obtained 
with the AD method, and confirms our previous result of the evidence of an important 
contribution of the isoscalar (7(500) meson in — > Tt~7t~^n^ decay [1]. We use the 
well known /2(1270) tensor meson in the crossing channel, as the base resonance, to 
extract the phase motion of the low mass kk scalar amplitude. We obtain a ^(^is) 
variation of about 180° consistent with a resonant a (500) contribution. We also obtain 
good agreement between the FSI jcorr observed with AD method and the 7 observed in 
the full Dalitz plot analysis. 
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